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Abstract 

We explore the bistability effect in a dimensionally scaled semiconductor nanostruncture con- 
sisting of a diluted magnetic semiconductor quantum dot (QD) and a reservoir of itinerant holes 
separated by a barrier. The bistability stems from the magnetic phase transition in the QD medi- 
ated by the changes in the hole population. Our calculation shows that when properly designed, 
the thermodynamic equilibrium of the scaled structure can be achieved at two different configu- 
rations; i.e., the one with the QD in a ferromagnetic state with a sufficient number of holes and 
the other with the depopulated QD in a paramagnetic state. Subsequently, the parameter window 
suitable for this bistability formation is discussed along with the the conditions for the maximum 
robustness/non- volatility. To examine the issue of scaling, an estimation of the bistabiity lifetime 
is made by considering the thermal fluctuation in the QD hole population via the spontaneous 
transitions. A numerical evaluation is carried out for a typical carrier-mediated magnetic semicon- 
ductor (e.g., GaMnAs) as well as for a hypothetical case of high Curie temperature for potential 
room temperature operation. 

PACS numbers: 72.20.Ht,85.60.Dw,42.65.Pc,78.66.-w 



The magnetism in semiconductors is the basis for the emerging field of spin-polarized 
electronics, or spintronics.— Substantial progress has been made during the past few years, 
particularly in the materials development. The advantages of the semiconductor-based sys- 
tems over the metallic counterparts include the controllability of the ferromagnetism (via the 
bias^ and/or doping) and the potential compatibility with the modern Si-based processing 
technology. 

Recently, a theoretical study further explored the opportunities offered by the electrically 
controlled magnetism. The calculation found that a properly designed structure consisting of 
magnetic and non-magnetic semiconductor quantum wells (QWs) can exhibit the bistability 
with respect to the paramagnetic-ferromagnetic (PM-FM) phase transition when the process 
is controlled by the itinerant carriers (holes). ^ This bistability effect in the two-dimensional 
system was predicted to persist even at temperatures nearly as high as the critical temper- 
ature Tc of the PM-FM phase transition. Subsequently, a non-volatile memory application 
was suggested citing the successful growth of transition metal doped semiconductors that 
are FM at or above room temperature (e.g., the nitrides,— Ge,-^ as well as some II-VI's^). 

For practical realization of the proposed device application, it is highly desirable to re- 
duce the size of the magnetic layer (i.e., the active part of the memory) without the loss of 
high temperature operability and non-volatility. Consequently, a magnetic semiconductor 
quantum dot (QD) that can exchange itinerant holes with a reservoir provides an interesting 
opportunity. However, unlike the QW case examined earlier,— the thermal fluctuation result- 
ing from spontaneous hopping between the two possible stable states cannot be neglected in 
the scaled structure due to the limited number of carriers populating the QD. Hence, it is 
the purpose of this paper to investigate theoretically the effect of the reduced dimensionality 
and explore the potential bistability conditions based on the electrically controlled magnetic 
phase transition in the magnetic semiconductor nanostructures. 

The speciflc structure under consideration consists of a single diluted magnetic semicon- 
ductor (DMS) QD separated from a large reservoir of itinerant holes, which controls the 
chemical potential /io of the system. For simplicity, we assume /io ^ ksT and ignore the 
possible temperature dependence of fiQ. A non- magnetic QW fllled with itinerant holes (for 
example, through the modulation doping, etc.) can be used as the desired reservoir (see 
Fig. 1). 

To accurately describe this system, one needs to know the energy structure of the QD 
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that is a function of multiple parameters. Particularly, the magnetic interactions must 
be taken into account that leads to the PM-FM phase transition in the DMS QD. Note 
that the analysis of different magnetic phase transition mechanisms is beyond the scope 
of the present study (see, for details, Ref. 0). Instead, we assume that the main magnetic 
properties [such as the critical temperature and its dependence on the hole concentration, 
and the magnetization dependence on temperature M = M{T)] can be obtained from the 
measurements of the relevant DMS. This allows the use of a semi-phenomenological approach 
to calculate the free energy of the system. 

In the present model, we approximate the free energy of the QD as the sum of two terms, 
the magnetic (Fm) and non-magnetic (-F/v) contributions. If the DMS QD is not far from 
the PM-FM transition, the Landau expansion over the magnetization M can be applied for 
Fm'- 



The parameters a, b and Tc are the functions of hole population j in the QD; particularly, the 
dependence Tc = Tc{j) plays a crucial role in the magnetic instability.^ In addition, a and h 
can be expressed in terms of the fundamental properties of the magnet: The Curie- Weiss law 
for magnetic susceptibility x = Co/(T — Tc) at T > Tc defines a = l/4Co, while spontaneous 
magnetization Mg = Mqa/I — T/Tc at T < Tc provides b = aTc/2M^. Since Eq. is 
assumed to fully describe the magnetization of the DMS QD, the thermodynamically stable 
state corresponds to the specific magnitude of M that gives rise to the free energy minimum 



where 6{x) is the Heaviside step function. It is important to note that the localized spin S of 
the magnetic ions provides the main contribution to the DMS magnetization, whereas that 
of the itinerant carriers add a minor role. The parameters Mq and Cq can be easily estimated 
for Nm localized spin moments leading to the estimation Mq/Cq = 3SNm/8{S + 1), which 
is independent of carriers. Thus, the only dependence of Fm on the hole population j comes 
from the term Tc = Tc{j) in our approximation. 

To obtain a numerical value for j, one must also incorporate the non-magnetic part Fn 
of the free energy for j particles located in the QD. Unfortunately, the calculation of Fn 
requires the approaches very specific to each individual case as it depends on the details such 
as the material composition, size and shape of the QD, presence of dopants and external 



a{Tc - T)M^ + bM^. 



(1) 





3 



fields, etc. Consequently, this problem cannot be solved in a general manner. To proceed 
further, we treat the QD as a scaled QW with a finite lateral size (and thickness); this is 
analogous to a nanodot embedded in a barrier. 

The potential profile of the sample structure along the growth (z) direction is schemati- 
cally illustrated in Fig. 1 from the hole representation, as is the case throughout the paper. 
It is convenient to split F/v into two parts, F/v = Ej + Fi(T,j). The first term, 

E,=jU + ^jij-l)C, (3) 

accounts for the energy acquired by j particles due to their localization in the QD with the 
ground state energy U as well as their Coulomb repulsion energy (C = e^/ey/Ao, where e is 
the electron charge, e the dielectric constant, and Aq the lateral cross-section of the QD).— 
Then, the remaining part, 

Fi(T,j)=f^(T,/ii)+j/ii(j), (4) 
is similar to the free electron gas contribution with the thermodynamic potential^i 

fi(T,/ii) = J]ln[l + e(^^-"")/^^^]. (5) 

n 

En and /ii represent the energy spectrum (with the quantum number n) and the chemical 
potential of the QD when the influence of the magnetic interaction is excluded (i.e., the 
non- magnetic version of the QD). Hereinafter, we consider the lateral dimension Ao of the 
QD to be relatively sizable so that the energy gaps in the discrete energy spectra are smaller 
than ksT. Then, the sum in Eq. © can be approximated by an integral with the density 
of states mAo/nh?. From the relation j = —dQ{T,fii)/dfii, /ii(j) is found to be 

fi,{j) = kBT\n{e^^ -1), (6) 

where ^ = nh'^ /mAokBT and m is the in-plane hole effective mass. Subsequently, the 
thermodynamic potential [Eq. (jS))] can be expressed in the form 

oo 

Q[T, /.i(j)] = 1 In [1 + (e^^- - l) e"^] dx. (7) 



Equations (jHl) and (jlj) along with /ii(j) and Q{T,fii) from Eqs. (jHl) and determine the 
non-magnetic part of the free energy, while the total free energy of the DMS QD is the sum 
F = Fm + -Fjv. 
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Now we take into account that the QD is in contact with a large reservoir providing 
two-way exchange of carriers through the potential barrier (see Fig. 1). This leads to the 
establishment of a unified chemical potential that coincides with /iq of the reservoir. Thus, 
the equation that determines the population of the QD takes the form 

f^Qoij) = /^o- (8) 

Note that figoij) 7^ f^iU) since both the non- magnetic and magnetic interactions are con- 
sidered for figoU)- Since the chemical potential figoij) of the QD can be expressed as 
f^Qoij) = dF/dj in general, the stable solutions of Eq. (jH)) must correspond to the local 
minima of F = F{j) or equivalently dfiQD{j)/dj > 0. 

Finally, the desired solutions require a specific expression for the dependence Tc = Tc{j) 
in Eq. (j2)). Following the experimental data of Ref. we adopt a semi-phenomenological 
description 

Te = r° (1 - e-"^«*) , (9) 

where is the asymptotic value of the critical temperature (at a sufficiently high hole 
concentration), t = T/T^, and a{= 1) is the fitting parameter that adjusts Eq. Q to the 
experiments. Subsequently, the QD population j can be obtained from Eqs. ©, ©, Q, 
and (IHI) as 



f^o-U = {j- 1/2)C + ksT In (e^^' - l) - 

T 
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For a numerical evaluation, let us assume the following set of parameters "typical" for a 
carrier mediated DMS (e.g., Gao.95Mno.05As): m = 0.13mo (mo is the free electron mass), 
e = 12.9, S = 5/2, Nm = 1.3 x 10^^ cm^3x(QD volume), and T° = 110 K. Figure 2 depicts 
the dependence [fJ^Qvij) — U] /ksT^ vs. j. Clearly, the results indicate that only one solution 
for j exists at sufficiently low or high energies U in reference to /io (e.g., dashed line 1 or 
3) corresponding to the only stable QD population at a given U . However, the moderate 
values of U can support multiple roots. In the case of dashed line 2, two of them (with 
the smallest and largest j) are stable considering the positive derivative {dfigo/dj > 0), 
while the intermediate solution is not {dfiQn/dj < 0). Note also that the stable solutions 
with the larger (smaller) j are realized in the FM (PM) phase of the DMS QD. Hence, this 
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demonstrates a bistable state for a properly designed QD in terms of the hole population or 
the magnetic phase. Figure 3 provides the ranges of U (assuming a fixed /iq) and T where 
the bistability can be expected in the system under consideration. Clearly, the calculation 
results indicates a large window in the system parameter space where the bistability is 
possible. The maximum operating temperature may be nearly as high as T°. 

To achieve a bistable state robust against the thermal fluctuations (i.e., non- volatile) 
for the possible memory application, it is necessary to select a condition that provides the 
maximal separation AF — min{Finax — -^Pmin, -^max — -^Fmin} between the local maximum 
Fmax and each of the local minima Fpmin and -FFmin (for the PM and FM phase of the DMS 
QD, respectively) of the free energy F{j). At a given temperature, one can find the QD 
potential alignment that results in the maximal AF. Figure 4 illustrates the behavior of 
F{j) at three different shifts A/i — iiq — U: curve 1 (3) represents the case of the monostable 
state in the FM (PM) state, whereas curve 2 exhibits two local minima in the FM and PM 
states, respectively, separated by the maximal AF. The temperature dependence of the 
maximal AF and the associated optimal potential shift A/x is plotted in Fig. 5 along with 
the mean values of the particle numbers jp, jm and jp corresponding to Fp min,Finax and 
Fp min- Our analysis shows that the "thermodynamic barrier" AF decreases drastically as T 
approaches [curve 2 in Fig. 5(a)]. Obviously, the system based on the PM-FM transition 
becomes much less stable near the critical temperature. 

Note that the mean value j is reached through the balance of the particle flux to and from 
the QD. Each of these incidents transfers one particle via the thermally activated processes 
with a characteristic time tq, which depends on the temperature, height and width of the 
energy barrier separating the QD and reservoir, etc. The lifetime is defined as the time it 
takes to develop a sufficiently large fluctuation to induce a transition from the state initially 
at Fpmin = F[jp) to that at Fpmin — F{jp)- Apparently if the system reaches the state 
at the local maximum Fmax = F{3m)-i further evolution can result in either the PM (with 
j — jp) or the FM (with j — jp) QD state with an approximately equal probability; hence. 
Tit can be estimated as the reciprocal probabihty for the process Fp 

min Fyiihx Or jp > jVrt- 

The formation of fluctuation Aj = jp — j,„t can be considered based on the sequential 
process of hole withdrawal from the QD. For the first hole transfer out of the QD, the 
characteristic time of this process is Tq as defined above. Due to the finiteness of the QD 
hole population, this reduces the chemical potential by A//(l) = /J'QdUf) — /J-qdUf — 1)- 
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Consequently, the new time constant becomes tq exp[A/i(l) / fc^T] when the next hole escapes 
from the QD, provided no particles are injected into the QD from the reservoir. Hence, 
the mean time necessary for the Aj fluctuation through the sequential withdrawal can be 
estimated as 

T^(Aj) = To 5^ exp[Afi{j)/kBT]. (11) 

j=0 

On the other hand, the probability of no hole injection from the reservoir during this time 
span Tu, is Pw = exp[—Tw{Aj)/To]. The frequency of appearance for such a rare occasion is 
Pw/tq and the desired lifetime is 

Tit = roexp[T^{Aj)/To]. (12) 

A similar expression is shown to apply to the fluctuations Fp-^^^ -^max or jp — > j^- 

One can see that the estimated lifetime (or the bit retention time) depends crucially on 
the operating temperature and the QD sizes, which determine the number of terms jm in 
the sum of Eq. (jllj) . Figure 6 analyzes the results for the GaMnAs QD with the dimension 
of 25 X 25 X 5 nm'^ and 15x15x5 nm'^ assuming = 110 K. As expected, the bistable 
state becomes less stable (i.e., shorter Tu) as the QD size shrinks. The thermal fluctuation 
clearly has a bigger impact in this case due to the finite number of holes in the QD. For 
the two structures considered (or those of similar sizes), a practically non- volatile condition 
(i.e., sufficiently long bit retention) may be achieved when operating below approx. 75 K. 
To = 1 ns is used for the calculation. 

To examine the feasibility of room temperature application, it is desirable to extend our 
consideration of magnetic semiconductors to those emerging ones with potentially much 
higher critical temperatures. Note for instance two recent reports of the DMS with > 
300 KM- Since the search for an ideal DMS is only at the beginning stage, we assume a 
hypothetical material with the characteristics similar to GaMnAs except in Eq. Q, which 
is treated as a variable in a wide temperature range. Figure 7 depicts the estimated lifetime 
as a function of T°, while T is fixed at 300 K. For a sufficiently long Tu in this case, the 
desired material needs of approx. 550 K or higher. 

In summary, we investigate theoretically the effect of the reduced dimensionality and 
explore the potential bistability conditions based on the electrically controlled magnetic 
phase transition in the magnetic semiconductor nanostructures. The analysis is based on 
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a semi-phenomenological model that assumes the common magnetic behavior and a simple 
hole energy spectrum in a DMS QD. When properly designed, the calculation predicts 
the possibility of controlled switching between the stable PM and FM states in the QD. 
The parameter window suitable for the bistability formation is discussed along with the 
conditions for the maximum robustncss/non- volatility. An estimation of the bistability 
lifetime limited by the thermal fluctuation provides a guideline for its potential application 
as a room temperature low-power, high-density memory element. 

This work was supported in part by the Defense Advanced Research Projects Agency 
and the SRC/MARCO Center on FENA. 
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FIG. 1: Schematic diagram of the structure containing a DMS QD (in two different phases) and 
a nonmagnetic quantum well (NQW) reservoir separated by a barrier. The energy orientation is 
provided from the hole point of view in the valence band. Left: The PM phase corresponds to 
disordered magnetic ion spins (small arrows) and the lack of the magnetic contribution to the 
hole energy. This is achieved when the holes (small circles) only weakly populate the QD with a 
discrete energy spectrum and, thus, cannot change its magnetic state. The ground state U of the 
QD is relatively high in reference to the chemical potential /io of the hole reservoir. Right: Another 
thermodynamically stable state (at the same external conditions) is possible when the magnetic 
ions are ferromagnetically ordered. Magnetic interactions can decrease the hole potential that the 
ground state of the DMS QD is now sufficiently below //q; the equilibrium hole population is high 
enough to stabilize the FM phase. Switching between the PM and FM states can be achieved by 
applying a gate bias Vg that populates or depopulates the DMS QD. 
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FIG. 2: Chemical potential of the DMS QD with the thickness of 5 nm and the cross-section of 
25 X 25 nm^ as a function of the QD hole population [Eq. (jlUj) ]. The parameters of Gao.95Mno.05As 
are assumed with T = 77 K as discussed in the text. The solutions of Eq. (|1U|) can be found as 
intersections of the solid curve with the horizontal line corresponding to a certain value of Afi 
(=/io — U). Two cases Afi/T^ =6 and 0.5 (dashed lines 1 and 3) provide monostable FM and PM 
states, while dashed line 2 (A/i/T^ = 2.7) depicts the bistable state. Stable solutions are indicated 
by single-head arrows and the unstable one by the horizontal double-head arrow. 
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FIG. 3: Phase diagram of the parameter space indicating the potential bistabihty region (the 
shaded area). The PM and FM denote the monostable areas corresponding to the PM and FM 
QD states, respectively. The same parameters as in Fig. 2 are assumed (T^ = 110 K). 
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FIG. 4: Free energy of the QD calculated as a function of hole population for three different values 
of A/x (=/xo — U): A/x/r^=6 (curve 1, dotted line); 2.7 (curve 2, solid line); 0.5 (curve 3, dashed 
line). The single minima of curves 1 (FM phase) and 3 (PM phase) correspond to the vicinities 
of the right and left boundaries of the bistable area in Fig. 3; curve 2 represent the bistable case 
with the optimal free energy barrier height separating two local minima. The same parameters as 
in Fig. 2 are assumed (T = 77 K, r° = 110 K). 
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FIG. 5: Maximal free energy barrier height (in units of ksT) between the local maximum and 
the closest local minimum in the bistable case (curve 2) and the corresponding optimal chemical 
potential (curve 1) as a functions of temperature. Curves 3—5 in the bottom pane shows the hole 
population at the FM minimum (curve 3), the PM minimum (curve 5), and the local maximum 
(curve 4) of F{j) for the bistable case shown in the top pane at the corresponding temperature. 
The same parameters as in Fig. 2 are assumed (T^ = 110 K) 
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FIG. 6: Bistability lifetime vs. temperature at T° = 110 K. Two different QD dimensions are 
considered: (1) 25 x 25 x 5 nm^ and (2) 15 x 15 x 5 nm^. The mean time tq of particle exchange 
between the QD and reservoir via thermal processes is assumed to be 1 ns. Other parameters are 
the same as in Fig. 2. 
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FIG. 7: Bistability lifetime vs. temperature for a hypothetical DMS material with the characteris- 
tics similar to GaMnAs except in Eq. Q, which is treated as a variable in a wide temperature 
range. T is fixed at 300 K and three different QD sizes are considered: (1) 25 x 25 x 5 nm^; (2) 
15 X 15 X 5 nm^; and (3) 25 x 25 X 3 nm . The mean time tq of particle exchange between the QD 
and reservoir via thermal processes is assumed to be 1 ns as in Fig. 6 
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